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Abstract 

The annihilation pd — ► 2tt~tt + p at rest is analyzed. Assuming that the 
deuteron wave function can be represented by the Fock's column consisting of 
different possible states ( NN, NN*, AA, NNir,...) some enhancement in the 
distribution over the invariant mass of the ir~Tr~ system in the mass region 
1.4 - 1.5(GeV/c 2 ) is predicted. This can be caused by the possible existence of 
a A — A component in deuteron. Using this prediction and new experimental 
data of the OBELIX collaboration (LEAR, CERN) one can estimate the upper 
limit for this exotic state of deuteron. 



PACS: 25.10.+S, 25.40.Ep, 24.10.Jv 

Keywords: Inelastic proton scattering, relativistic models 



2 



Over the past decades the deuteron structure at short distances has been investi- 
gated rather intensively. There are many theoretical approaches trying to study this 
problem. For example, within the quark model picture the admixture of 6q-components 
has been suggested §]. This corresponds to an admixture of baryon resonances 
within the asymptotic basis of baryonic states (JV(938), A(1232), JV(1440), ...) §, ||. 
New information about the existence of, for example, iVW* (1440) in a deuteron can be 
obtained by studying subthreshold antiproton production in pd and dd reactions 0. It 
was shown that the probability of this component in a deuteron can be less than 0.5%. 
Attempts to detect baryon resonances in a deuteron were made by many experimental 
groups H |7j and || |[|. As a rule, pion, photon, neutrino and proton beams of interme- 
diate energy together with fixed target setups were used for this purpose. These groups 
have tried to demonstrate the presence of N* in nuclei by shaking loose an iV* that 
preexists in the nucleus acting as a spectator in some reactions. The main problem 
for such experiments was to distinguish between the signal from the internal N* and 
the signal from another states produced by the conventional way, e.g. via rescattering. 
The experimental results are given as upper limits, ranging in (0.4 — 0.9)%. 

We propose alternatively to explore antiproton annihilation on the deuteron to 
provide new information on the possible existence of the above mentioned states in 
the deuteron. Here we are intend well known measurements of pd annihilation at rest 
performed by the OBELIX collaboration (see e.g. jUj). New data obtained by this 
experiment |nj show the very interesting behavior of the 7r~7r~ system invariant mass 
distribution in reaction — > 2tt~7t + p. In this paper we concentrate on detailed analysis 
of this annihilation channel. We show also that some enhancement in the spectrum 
over the invariant mass of two 7r~-mesons observed by OBELIX can be explained by 
the possible existence of a A — A component in the deuteron. 

The general expression for the differential cross section dafdm^-K for the above 
mentioned annihilation reaction in flight within the framework of the impulse approx- 
imation, see Fig.l, can be written in the following form: 

d(J < 2tTYli i 2'n- f j 7-1 / 2 7i /r1\ t~i / 2 2\ /i \ 

dm 27T (27r) 2A 1 / 2 (s, m 2 ,Mj) J 

where the following notation is introduced: is the invariant mass of two pions; 

A 1 / 2 (x, y,z) = (x — (y/y — \fz) 2 ) l l 2 (x — (v/y+ ^/z) 2 ) 1 ^ 2 ] s is the square initial energy in 
the p — d c.m.s.; is the square invariant mass of three pions; m, are the masses 
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of the nucleoli (antinucleon) and the deuteron, respectively; m n is the mass of the 
neutron inside the deuteron which, in principle, is off-shell; R 2 (s,m 2 ,M^) is the two- 
particle relativistic invariant phase-space for the binary process p + d —>■ 3tt+p when in 
the final state there are a proton and a particle with the invariant mass of three pions; 
R^s^, w? nl m 2 ) is the three-particle relativistic invariant phase-space for annihilation of 
the initial antiproton on the intra-deuteron neutron pn — > 27r~7r + ; Tp d ^ 37rp (s, s 37r ,t,ti) 
is the amplitude of the process pd — > 2n^n + p, t = (pd — P P ) 2 , where Pd,P P are the 
four-momenta of the initial deuteron and final proton, respectively; t\ = (p p — p n ) 2 , 
where pp, p^ are the four-momenta of the p and of the final pions, respectively. 

Within the framework of the impulse approximation, see Fig.l, this amplitude can 
be written in the factorized form: 

Tpd->3n P {s, S 3n , S 2 -7r,t, *l) = 0d(Pp( S ) *) )fpn^2ir~TT+ (S3n, S 2n , h) (2) 

where (pd and fp n ^2n-ir+ are the deuteron wave function and the amplitude of the 
annihilation process pn — > 27r _ 7r + ; p p is the three-momentum of the final proton. Now, 
inserting now expressions for R 2 , R 3 |L2| and Tp d ^ 3np given by eq.(§) into eq.(l) one can 



get the following form for the distribution over the invariant mass of the 7i + n -pair: 
da 2m 7T + 7T -n 3 I f , f dt 



dm^- 2(2tt)«2A( S , m 2 ,Mj) 2W 2 {s^- ^ 2 ^ 2 ) J **** J 2A 1 /2( Ss ^ m2 ) ( 3 ) 

I </>d(p( S (0) f J dtl J ds 27T -\ fpn->27T-Tr+(s 37T ,S 27T ,t 1 ) \ 2 

A similar expression can be written for the distribution over the invariant mass of the 
7r~7r~-pair: 

da 1 2m 27T -7Y 3 1 f , f dt . . 

ds ^ / 0M/9/ 1 5T ( 4 ) 



dm^-v- 2(2n) s 2\(s,m 2 ,M 2 )2\ 1 / 2 (s 27T -,fi 2 ,n 2 ) J J 2\ 1 / 2 (s 37T , t, m 2 , 

I f J dt l J ds w+7r - \ /p„_ 2 7r-7r+(S37r, S 2 ^,tl) f 

where s 2n = m\ % is the square invariant mass of 2ir~, fi is the pion mass. For calculation 
of the spectra given by eqs. (3), (4) the amplitude of the annihilation pn — > 2h~'k + 
should be known. 

Actually, the question arises how to calculate the amplitude fp n ^2-K--K+- It was 
parametrized by a form that resulted in a satisfactory description of the experimental 
data for dN/ dm^^- . According to the OBELIX data [fl ij the spectrum dN/ dm n + n - 
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fpn^2n-TT+(s 37 T, S 27T ,ti) \ 2 = B 2 (s iv + lv -) — 1 ^ F 2 (s 37r ) n - j ■ _ 2 \2 ^^' 1 ) (^) 









(1*1 


+m 2 ) 2 



has a resonance-like shape, see Fig.2, corresponding to the possible production of p- 
meson at m w + w - ~ 0.77(GeV/c 2 ) and / 2 - m eson at m n + n - = 1.26(GeV/c 2 ). The de- 
pendence of | fp n ^2n-TT+( s 3n, S2n,ti) | 2 on Sp n = s 3n can be found using the model 
of the Reggeized meson exchange in the t-channel of the reaction 7i + n — ► 7t + tt~p, 
where t > 0. and t = (p„ — p p ) 2 = (p„ + pp) 2 = Sp n = s 3w . The ^-dependence of 
I fpn^2n-Tr+{ s 3Tr, S2td*i) | 2 can De found using the usual one-baryon exchange model of 
the reaction pn — > 2n^n + . Therefore, the square of the amplitude fp n ^2-n-n+(s3n, s 2 -k, t\) 
has been parametrized by the following form: 

[s 3n + p 2 

where -Fi(s37r) is the pion form-factor corresponding to the NNn- vertex by the Reggeized 
meson-exchange of ir + n — > 7T + 7T~p process, it can be parametrized by the usual form 
Fi(s 3v ) = exp{—s 3n R\) where the parameter R 2 ~ 2GeV~ 2 |13| is used; F 2 (ti) is the 
baryon form-factor corresponding to the iViW-vertex by the one-baryon exchange in s- 
channel of the reaction pn — > 27r~7r + , it was taken in the form F 2 (t 1 ) = A 2 B /(\ t\ | +A. 2 3 ) 
where A# = 0.4 — 0.5(GeV/c) 2 fll4| ; B(s n + n -) is the function corresponding to the res- 
onance and nonresonance mechanisms of the annihilation process pn — > 27r~7r + , it was 
parametrized by the following form: 

B{s^-) = Wp — -- + Wh — -- + iw (6) 

y/Sn+w- -m p - iTp/2 ^/s^+~ - m h - iT h /2 

Really, the form (^j) is some parametrization of the square of the amplitude of the 
process pn — > 2ir~7T + , and application of the Reggeized one-boson exchange model can 
be considered only as a prompt to find the parametrization form (|5|). 

Now, inserting eq.(B|) into eq.(3) one can calculate the spectrum da / dm n +^- . Ex- 



perimental measurements of this spectrum were performed at rest and presented in [ |TT 
Therefore, the spectrum given by eq. (3) was calculated at a very small initial antiproton 
energy satisfying to the experimental resolution and normalized by the total number of 
events, N, corresponding to experimental measurements. This spectrum dN/dm 7r + n - 
and the experimental data are presented in Fig.2, the parameters w p ,Wf 2 ,w and R 2 
were found by fitting the OBELIX experimental data JO] for dN/dm^+T^-, they are 
w p = 1.2, w h = 1.35, w = 35.0 and R\ ~ 2.(FeV/c)~ 2 . 

Now, inserting the same form for | fp n ^2n~n+{s3TT, s 2lT , ti) | 2 given by eq.@ into 
eq.(4) one can calculate the spectrum dN/dm 2n - which is presented in Fig.3 (the 



dashed curve) together with the experimental data |TT . 
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J ds n + p R 2 (s, s 27T -, s w+p )R 2 (s 2n - , /i 2 , /i 2 ) (7) 



Now, consider the problem related to the possible existence of baryon resonances 
in a deuteron. There is a point of view |15|, [16j that a deuteron wave function (d.w.f.) 
can be considered a Fock column the lines of which are NN, AA, NN* , NNw, etc., 
components. There are theoretical models within the framework of which the d.w.f. is 
considered taking into account these baryon resonances in a deuteron. We used one 
of these models [|l], |2j in order to estimate the upper limit of the probability for a 
AA-component to exist in the deuteron. The possible mechanism of the annihilation 
pd — > 37rp is depicted in Fig. 4. 

The contribution to the spectrum da/dm 27T - corresponding to this graph of Fig. 4 
can be presented as: 

da AA _ 2m 27r - 
dm 2n - ~ (27r) 8 2A 1 /2( S;m 2 ;M 2 

^(-Stt+p, W 2 , /i 2 ) | Tpd^Tr-n+p^S, S 2n - , S n + p , ti) | 

The amplitude Tp d ^ 27T ~ n + p (s 27T - , s w +p, t\) corresponding to the graph of Fig. 4 is written 
in the following form: 

where 4> aa (pa) is the part of the d.w.f. corresponding to the A A component of Fock's 
column; fpA-^2w- is the amplitude of the process pA~ — > 2n~] fA++^TT+ P ( s -K+p) is the 
function corresponding to the down vertex of Fig. 4, e.g., the decay of A ++ — > 7i + p, and 
depending on the square of the invariant mass of the 7t + p system s^+p ■ The form of 
I fpA-^2n- 1 2 can be found using the one-baryon exchange model, e.g.: 

I fpA~^2n- | 2 (tl) = / 1 . | ] o\2 F 2 ( 9 ) 

(| ti | +m 2 ) 2 

The function fA++^TT+p( s n+ P ) has the Breit-Wigner form and can be written as the 
following: 

C 

I l 2 = (v ^_ mA)2 + ri/ 4 (10) 

where tua is the A-isobar mass, Ta is its width, C is some constant. The form of the 
d.w.f. 4> aa (pa) taking into account the possible existence of a AA-component in the 
deuteron was taken from |l|, 0] obtained within the framework of the quark model and 
had the following form: 

Ah 2 

<P AA (Pa) = v / ^( 0F(3/2)3/2 ) 1/2 ^(-^/3) (11) 
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where sa is the so-called spectroscopic factor 0; b is the slope parameter 



Now, inserting R 2 (s, s 2n - , s n + p ), R 2 (s 2n - , fi 2 , ^ 2 ), R 2 (s n + P ,m 2 , fx 2 ) and eqs.(|8|-[lTD into 

rAA , 



eq.(7) one can get the expression for the contribution da j dm 27T - : 

da AA 2m 27r -n 3 r X 1/2 (s 27T -, m 2 , n 2 ) 

dm 2w - (27r)82A( S ,m2,Mj)2A 1 /2( S27r _^2 5/i 2 ) J s - + p 2s , +p 1 ' 

Then, the total spectrum da/dm 2w - within the framework of the impulse approxi- 
mation can be obtained as the incoherent sum of da^ / dm 27T - and da AA /dm 2n -, e.g.: 

da da 1 | da AA 



dm 2lT - dm 2n - dm 27r 
The corresponding spectrum dN / dm 27T - calculated in the same manner as dN / dm^n 



and the OBELIX experimental data [11] are presented in Fig. 3 (solid line). By the 



calculation of da AA / dm 2n - the joint parameter value sa * C 2 has been found fitting 



the experimental data for dN/dm 27T - | IT| . According to our calculations the spectrum 
da AA jdm 2l ,- may contribute to the total spectrum at m 27T - = 1.4 — 1.5(GeV/c 2 ). 
Therefore, the small enhancement in dN/dm 27T - observed experimentally at m 2n - = 
1.4 — 1.5(GeV/c 2 ), see Fig. 3, can be interpreted as the contribution of da AA / dm 2n - . 

At this point the question arises about the contribution of corrections to the im- 
pulse approximation of the annihilation process pd — > 2tt~tt + p, see Fig.l. As for 
this possible correction, the two-step mechanism [|17]] can contribute to the spec- 



trum da/dm 2lT . Our estimations and the results of ref.[17] showed that this contri- 
bution is very small for the invariant mass value of 7r + p m^+ p < 1.1 — 1.15(GeV/c 2 ) 
tha corresponds to m 2n - > 1.3 — 1.35(GeV/c 2 ) and can't give any enhancement at 
m 2n - = 1.4 — 1.5(GeV/c 2 ) to the spectrum dN/dm 27T - , when according to |lTJ the mo- 



mentum of the proton p > 400 MeV/c. In principle, the two-step mechanism can give 
some sizeable contribution to this spectrum at m 2n - ~ 0.8 — l.l(GeV/c 2 ). However, it 
is very difficult to calculate it accurately because there is a large uncertainty related 
to the off-shellness of the virtual meson in the intermediate state 0, |i~9| . Thus, ac- 
cording to the above this AA component can contribute in the high mass kinematical 
region where the two-step mechanism of the annihilation process pd — > 2tt~7t + p can be 
neglected. 

So, let us assume that a small enrichment in the n~ tt~ system invariant mass 
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distribution in the reaction pd — > 27r~7r + p is a manifistation of the AA component of 
the deuteron. 

In this case, due to the lack of statistics, we can only obtain an upper limit for 
the process of Fig.4. Fig.5(a) shows the scatter plot for the invariant mass of the 
7r~7r - system versus the invariant mass of the n + p for pd — ► 27r _ 7r + p annihilation via 
the channel of Fig.4. There is a some small enrichment in the upper part of the n + p 
system in the A ++ (1232) region. In Fig. 5(b), the invariant mass distribution of the 
7T~7r~ system in this reaction is shown. The solid line corresponds to all the selected 
events. The hatched histogram corresponds to the events with the invariant mass of 
the n + p system satisfying the following condition: (M^+p — M A ++\ < 120 Me V. As 
one can see, there is a small bump in the high mass part of the 7r~7r~ invariant mass 
distribution, but this bump has a small statistical significance, so we can only estimate 
the upper limit for the reaction of Fig.4. To do this the distribution in Fig. 5(b) was 
approximated at a first stage by a second order polynomial function. The parameters 
of this polynomial function were fixed and the experimental invariant mass distribution 
of the 7r~7r~ system was fitted at the second stage by the function: 

3 

F = A 1 x BW(T v - v -,M v - v -,m) + A 2 x ^ b, x to*" 1 (14) 

i 

where BWiT^-^-^ M n - n ~, m) is the Breit-Wigner function. During the fit, the param- 
eters A\, A2, Mtt-tt- were free, the parameters of the polynomial function were fixed 
from the first stage of approximation and the T^-^- parameter was fixed at a value 
of 60 MeV according to theoretical prediction. Using the results of this fit, the upper 
limit on the reaction of Fig.4 was estimated to be: 

Yp (A-A+ + )^2n-n+p < 6.5 X 10~ 5 (15) 

with a 90% confidence level. 

Finally, let us suppose that the branching ratio of the reaction pA~ — > 2n~ is 
approximately the same as for pN — > 2tt, that is ~ 5. x 10~ 3 Using this 
assumption and the upper limit we can roughly estimate the probability to find 
the deuteron as the AA configuration. This probability is 

Paa < 1-% 

Note, that in the near future a factor of 10 higher statistics will be available for analysis 
from the OBELIX collaboration. 
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Let us now make a short conclusion on our investigation. Some enrichment of the 
scatter plot for m 27r - versus m^+p for the pd — > 2ir^7Y + p process and the corresponding 
small enhancement in the spectrum dN/drri2 7T ~ at m2 V - = 1.4 — 1.5(GeV/c 2 ) observed 
experimentally |Tl| can be caused by the possible existence of a A — A component in 
the deuteron. The analysis of this phenomenon shows that it is difficult to calculate 
the absolute value of its probability. But the theoretical prediction about the width 
value of this enhancement allows us to estimate the upper limit for the A — A exotic 
state of the deuteron. 

The authors wish to thank V.G.Neudatchin, I.T.Obuchovsky, I.Chuvilsky and also 
M.P.Bussa for valuable hints and inspiring discussions. 
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